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ABSTRACT 

We present N 2 D+ 3-2 (IRAM) and H 2 D+ In - 1 10 and N 2 H+ 4-3 (JCMT) maps of 
the small cluster-forming Ophiuchus B2 Core in the nearby Ophiuchus molecular cloud. 
In conjunction with previously published N 2 H + 1-0 observations, the N 2 D + data reveal 
the deuterium fractionation in the high density gas across Oph B2. The average deu- 
terium fractionation Rd = A r (N 2 D + )/A^(N 2 H + ) ~ 0.03 over Oph B2, with several small 
scale Ro peaks and a maximum Rd = 0.1. The mean Rd is consistent with previous 
results in isolated starless and protostellar cores. The column density distributions of 
both H 2 D + and N 2 D + show no correlation with total H 2 column density. We find, 
however, an anticorrelation in deuterium fractionation with proximity to the embedded 
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protostars in Oph B2 to distances > 0.04 pc. Destruction mechanisms for deuterated 
molecules require gas temperatures greater than those previously determined through 
NH3 observations of Oph B2 to proceed. We present temperatures calculated for the 
dense core gas through the equating of non-thermal line widths for molecules (i.e., N2D + 
and H^D" 1 ") expected to trace the same core regions, but the observed complex line struc- 
tures in B2 preclude finding a reasonable result in many locations. This method may, 
however, work well in isolated cores with less complicated velocity structures. Finally, 
we use Rd and the H2D + column density across Oph B2 to set a lower limit on the 
ionization fraction across the core, finding a mean x e n m > few x 10~ 8 . Our results 
show that care must be taken when using deuterated species as a probe of the physical 
conditions of dense gas in star-forming regions. 

Subject headings: ISM: molecules - stars: formation - ISM: kinematics and dynamics - 
ISM: structure - radio lines: ISM 



Introduction 



The Ophiuchus molecular clou d is a source of ongoin g low-mass clustered star formation rela- 
tively ck2se_to_the Siinfce ~ 120 pc: lLombardi et al.ll2008l ). Through low resolution DCO + observa- 
tions, lLoren et al.l (|1990l ) discovered a number of dense Cores (labelled A through F) in the central 
Ophiuchus region with masses of a few to several tens of solar masses, which represented a range in 
evolutionary status given the presence or absence of embedded protostarfl Large-scale submillime- 
tre and millimetre observations of Ophiuchus identified several new Cores, and reveal ed the Cores 



were further fragmented in t o mu l tiple dense clumps , each with masses M < 1 M Q (jMotte et al 



19981 : 1.Tohnstone et al.l 12000 . 12004 : lYoung et al.l l200fil : ISadavov et al.1 l2010l ) . Analysis of the clump 



kinem atics suggests most are in virial equilibrium and hence not transient objects ([Andre et al 



20071 ). Recent Spitzer Space Telescope studies of infrared emiss ion towards the Cores have charac- 



terize d in detail the starless and protostellar clump population (jJorgensen et al.ll2008l : lEnoch et al 



2009). These objects are thus ideal locations to study the physical characteristics of dense cores in 
the process of forming small stellar clusters. 

To characterize the physical properties of the dense, cluster forming Ophiuchus Cores, we 
have previously stud ied Oph B, C and F through high resolution NH3 (1,1) and (2,2) observations 
(jFriesen et al.ll2009l . hereafter Paper I), and Oph B in both single-dish and interferometer N2H + 



1 In this paper, we describe Oph A, B, C, etc ., as 'C ores' since this is how these features were named in DCO + 
observations of the L1688 region bv lLoren etaD (|l990h . Since then, higher-resolution data such as those described 
in this paper have revealed substructure in these features that could be themselves precursors to stars, i.e., cores. 
To avoid confusion, we refer the larger features in Oph as 'Cores' and Core substructure identified by CLUMPFIND as 
'clumps'. 
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1-0 emission (jFriesen et al.l 120101 . hereafter Paper II). The NH3 and N2H 4 " observations probe gas 



at densities n ~ 10 3 4 cm 3 and n ~ 10 5 cm 3 , respectively. 



In Paper I, we identified (through the 3D clumpfind structure-finding algorithm; IWilliams et al 



1994 ) a number of individual clumps in NH3 emission in the Cores, which did not correlate well 
with those identified in continuum emission. We determined that Oph B and F were warmer on 
average than typically found for isolated cores (mean T% ~ 15 K), and we found no trend with tem- 
peratures and H2 column density calculated from submillimetre continuum emission. In contrast, 
Oph C is colder, with temperatures decreasing to a minimum of 9 K at the dust continuum emission 
peak. Line widths in both the Oph Bl and B2 sub-Cores were found to be wide and transonic, de- 
spite the presence of three embedded protostars in B2 and none in Bl. Oph F, with four embedded 
protostars, has regions of small and large line widths and complicated line structures. Again similar 
to isolated objects, NH3 line widths in Oph C are extremely narrow and consistent with being due 
to purely thermal motions. In an interesting result, the fractional abundance of NH3 with respect 
to H2, X(NH3), appeared to decrease with increasing H2 column density, A^FLi), suggesting NH3 
may be depleted towards locations of high column density and thus may not be tracing the densest 
core gas. Finally, Paper I revealed that cores with characteristics of both isolated and clustered 
environments can coexist in a clustered star forming environment. 

In Paper II, we again used clumpfind to identify small-scale structure in N2H + emission 
from Oph B and found that the N2H + clumps matched well NH3 clumps from Paper I, and were 
consequently again offset from continuum objects. In both studies, little difference was found in 
the gas properties towards embedded protostars compared with the general core gas. Line widths 
in Oph B2 remain large at the densities traced by N 2 H + . In Bl, however, N2H + line widths are 
significantly narrower than found in NH3 emission. Applying the gas temperatures determined in 
Paper I, we found that motions in Oph Bl were subsonic with a mean ratio of non-thermal line 
width to sound speed a NT /c s = 0.7 across the core, while in B2 a NT /c s = 1.3. A clear trend of 
decreasing A(N 2 H + ) with increasing iV(H 2 ) was found towards Oph B2, with emission suggestive 
of a N 2 H + hole towards the continuum emission peak found in high resolution N2H + 1-0 emission 
observed with the Australia Telescope Compact Array. 

The results of NH3 and N 2 H + observations towards Oph B2 in Paper I and Paper II indicate 
that the highest density gas may best be probed by molecular emission lines excited at higher 
critical densities, or by deuterated molecules which are expected to be abundant in cold cores 
where significant depletion has occurred. In cold molecular cloud cores, the deuteration process is 
initiated by the reaction 



H+ + HD ^ H 2 D+ + H 2 + AE 



(1) 



which is exothermic in the forward direction with AE = 232 K ([Miliar et al.lll989l ). At very low 
temperatures (< 20 K), the forward reaction dominates, making H2D + a key molecular ion in the 
enhancement of deuterated species in a molecular core. CO, which can destroy both H2D + and 
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Hg" via proton-transfer reactions, is additionally expected to deplete from the gas phase through 
freeze-out onto dust grains. The deuterium fractionation can thus be elevated orders of magnitude 
above the cosmic [D]/[H] ratio of 10 -5 for H3", as well as for other species which react with its 
deuterated forms. 

The ortho-H2D + In — 1 10 transition at 372 GHz was first d etected towards the low mass young 
stellar object (YSO) NGC 1333 IRAS 4A bv lstark et all jl99gh . Multiple detections of H?D+ and 



doub l y-deuterated Hj~ , D2H" 1 ", have follow ed in both starless and protostellar cores (jCaselli et al 



2003; Stark et al. 2004 



Vastel et alJbotMh. ey of H2D + emission towards 

sixteen objects by ICaselli et al.l (120081') . IWalmsley et al.1 (|2004l ) show the H2D + abundance de- 



termined for the dark core L1544 by ICaselli et al.l 120031 is consistent with a 'complete depletion' 
model, where heavy elements such as C, N and O (and, consequently, molecules containing these 
elements) are missing from the gas phase at the highest core densities. In this scenario where no 
heavy elements are left in the gas phase, H 2 D + and multiply-deuterated forms of are the only 
gas tracers of the physical conditions of the dense core. 

Based on our previous work in Paper I and Paper II, Oph B2 is a good target to investigate 
the distribution of deuterated species in a more complex star-forming core. To this end, we present 
in this paper N2D + 3-2, N2H + 4-3 and H2D + In — lio emission maps of Oph B2. In the following 
sections, we present the observations in ^21 and discuss the observed distributions of N2D + 3-2, 
H2D + In — lio an d N2H + 4-3 in $31 We fit the observed spectra with a single Gaussian (for H^D" 1 ") 
or a multiple-component Gaussian (for N2D + and N2H + due to their hyperfine line structure) and 
discuss the line velocity centroids and widths in §U and also present column density and fractional 
abundance calculations. In £}5] we look at the non-thermal line widths in Oph B2 as a function of 
critical density of the observed emission lines, and discuss trends in the deuterium fractionation. 
Finally, we present a lower limit to the electron abundance in B2, and summarize our results in SJBJ 



2. Observations 

In the following section we discuss the observations presented in this paper. We list the species, 
transitions, and rest frequencies observed in Table [H 



2.1. N 2 D+ at IRAM 

An On-The-Fly (OTF) map of the N 2 D+ 3-2 line at 231.321GHz in Oph B2 was made at 
the Institut de Radio Astronomie Millimetrique (IRAM) 30 m Te lescope during the 2 007 and 2008 
winter semesters using the HEterodyne Receiver Array (HERA, ISchuster et al.ll2004l ). HERA is a 
multi-beam spectral line receiver with nine dual polarization pixels arranged in a 3 x 3 array, with 
pixels separated by 24" [approximately twice the 231 GHz 11" beam full width half power (FWHP)]. 
The maps were performed with the array tracking the sky rotation, and rotated by 9.5° with respect 
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to the equatorial system, such that the spacing between OTF lines was 4". Observations were 
frequency-switched in-band with a throw of 2 x 6.9 MHz. The Versatile Spectrometer Assembly 
(VESPA) autocorrelator was used as the backend, with a 20 MHz bandwidth and 40 kHz channel 
spacing, corresponding to 0.05 km s -1 per channel at 231 GHz. 

The data were calibrated to T* A units at the telescope using the Multichannel Imaging and 
Calibration Software for Receiver Arrays (MIRA). Further calibration, including baselining and 
folding of the frequency-switched spectra, was performed using the CLASS software package!- The 
data were smoothed by two spectral channels along the frequency axis to improve sensitivity, for 
a final spectral resolution of 80 kHz, or 0.1 km s -1 . A second-order polynomial baseline was fitted 
to the non-line channels and subtracted from the individual spectra before folding the data. The 
forward (F e ff) and beam (B e ff) efficiencies of the telescope were interpolated from determined 
values at 210 GHz and 260 GHz, giving F eff = 0.92 and B eff = 0.58 at 231 GHz. The data were 
then calibrated to main beam temperature via the relation Tmb = {Feff I Beff) T\. For comparison 
with the N 2 H + 1-0 results in Paper II, we further converted the data to units of Jy /beam -1 and 
convolved the map to a final angular resolution of 18" to match the N 2 H + data. 

As the observations were performed in a pooled mode, the 4' x 3' map was observed in three 
4' x 1' strips to ensure good sensitivity in potentially limited time. Two-thirds of the map was 
completed in this manner in winter 2007, while the remaining strip was observed in winter 2008. 
The sensitivity of the map is consequently not entirely uniform, with the rms noise in the most 
southern 4' x 1' strip greater [0.15 K (Tmb)] than in the rest of the map [0.08 K (Tmb)]- Most 
emission was found in the central and northern map sections, however, so this noise increase does 
not affect our results significantly. 



2.2. H 2 D+ and N 2 H+ at the JCMT 

Simultaneous observations of the H 2 D + lio — In and N 2 H + 4-3 emission lines at 372.421 GHz 
and 372.673 GHz, respectively, were performed at the James Clerk Maxwell Telescope (JCMT) over 
the 07A - 8B semesters using the 16-receptor Heterodyne Array Receiver Program B-band receiver 



(HARP-B, ISmith et al.ll2008l ). HARP-B is a 4 x 4 pixel array, with array spacing of 30" and a 2' x 
2' footprint. The JCMT beam at 372 GHz is ~ 13" FWHM. We used the ACSIS correlator as the 
backend, configured to have a 400 MHz effective bandwidth (e.g., wide enough to encompass both 
lines), tuned halfway in frequency between the two lines, and a frequency resolution of 61kHz, 
or 0.05 km s -1 at 372.4 GHz. Observations of five separate HARP footprints were performed in 
position-switching mode to create an undersampled map of H 2 D + and N 2 H + over Oph B2. The 
footprints were placed such that the final beam spacing on the sky is 21" with a coverage of ~ 4' 
x 2'. The fifth footprint provided additional coverage towards the location of peak continuum 



2 see |http: / / www.iram.fr/IRAMFR/GILDASI 
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emission in B2, but has lower sensitivity compared with the other footprints. The observations 
were performed in band 1 weather only (T225GHZ < 0.05) and generally while Ophiuchus was above 
30° in elevation. The mean rms across the H2D + map is 0.04 K (TJ) and 0.05 K (TV) across the 
N2H + map. The main beam efficiency factors, tjmb, of the HARP receptors at 372 GHz are not 
currently well-known. At 345 GHz, the mean ijmb = 0.60, with an rms variation of ~ 5 % between 



detectors ([Buckle et al.ll2009l ). We therefore estimate t\mb ~ 0.6 with an uncertainty of ~ 10%. 



Data reduction was performed using the Starlink software package. Each night's data were 
first checked for noisy or malfunctioning pixels, and these pixels were flagged individually for each 
integration. A linear baseline was then removed from each integration, and the data were combined 
into a final cube with 15" pixels. To increase sensitivity, the data were averaged along the velocity 
axis by 2 channels, giving a final velocity resolution of 0.1 kms -1 . 



3. Results 



We first discuss the integrated intensity of N2D + 3-2, H2D + In — lio and N2H + 4-3 over Oph 
B2 and compare the species distributions with submillimetre continuum emission and the locations 
of embedded protostars. 

In Figure [Hi, we show the integrated N2D + 3-2 intensity (convolved to 18" resolution as de- 
scribed in §2.ip in Oph B2 against contour s of 850 continuum emission observed at the JCMT 
at 15" FWHM resolution (original m ap by I Johnstone et al.ll2000i , combined with a ll other SCUBA 
archive data by lJorgensen et al.l l2008. following; the description inlKirk et al.ll2006r ). Also identified 
are the locations of 850 ^m continu um clumps ( Jeirgensen et al. 20081 ) and Class I protostars [clas- 



sifications and locations taken from Enoch et al 



20091 . two of which were previously identified as 



Elias 32 and 33 (E32 and E33) bv lEIiaell 19781 : also VSSG 17 and 18 bv lVrba et alJll973 |. The N 2 D+ 
emission is confined to the northern edge of B2, and generally follows the continuum contours in 
this region, but a very close correspondence is only seen towards the north-east tip of B2. Three 
maxima of integrated intensity are present, of which two are co-located with continuum clumps, 
while the third is seen to the north-west of the continuum peak. The N2D + emission is weak 
towards the continuum emission peak, and avoids entirely areas near the protostars. 

In Figure [lb, we plot the N2D + 3-2 integrated intensity contours over a map of N2H 4 " 1-0 
integrated intensity, observed with the Nobeyama 45 m Telescope at 18" resolution. The N2H + 
observations and analysis were presented in Paper II. We also show the locations of N2H + clumps, 
identified through CLUMPFIND in Paper II. Overall, N2D + emission is found within the extent 
of N2H + 1-0 emission. In a similar fashion to the continuum emission, however, we find offsets 
between the locations of peak N2D + and N2H + integrated intensity. In particular, the N2D + 
integrated intensity maximum is offset from the integrated N2H + 1-0 maximum in B2 by ~ 20" to 
the northeast, ~ 1 beam width. The N2H + 1-0 integrated intensity maximum is itself offset to the 
east from the continuum peak by ~ 22" (Paper II). The northeastern N2D + integrated intensity 
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peak, which corresponds well with continuum contours, is also offset from the integrated N2H + 
intensity peak by ~ 20" to the east. 

We show in Figure [2] H2D + In — Iiq spectra across Oph B2 with overlaid 850 /xm continuum 
contours. We find that significant H^D" 1 " emission is present over much of Op h B2, with a larger 



extent than found in H2D" 1 " observatio ns towards any other cores so far observed (jCaselli et al.ll2003l : 



Vastel et al.l 120061 : iPagani et al.l 12009) . The H2D + emission again follows generally the continuum 
emission but avoids the embedded protostars in B2, similar to N2D + . While relatively strong 
emission is seen towards the continuum peak, the H2D + integrated intensity maximum (not shown) 
is also offset the continuum peak by an estimated ~ 20 — 30" (due to the map's undersampling we 
are unable to determine exactly the location of maximum integrated intensity). In contrast to the 
N2D + emission, the H2D + integrated intensity maximum is found to the northwest rather than 
the northeast (although N2D + emission is also found towards the H2D" 1 " peak). The pixel-to-pixel 
variations seen in Oph B2 may not be significant, however, due to the ~ 5% variation deter mined 



for the main beam efficiency factors between different HARP receptors (IBuckle et al.ll2009l ). We 
detect some E^D" 1 " emission in north-east B2 where was found strong N2D + emission, but do not 
see a significant integrated intensity maximum. (Note, however, that the HARP footprints do not 
extend past an R.A. of 16:27:35, and are therefore not sensitive to any emission at the tip of the 
continuum contours in this region.) Overall, where N2D + emission appears strongest in the north 
and east of Oph B2, H2D + emission is strongest towards B2's western edge. 

In Figure El we show N2H + 4-3 spectra across Oph B2 with overlaid 850 /jm continuum 
emission contours. We find bright N2H + 4-3 emission over most of B2. The greatest integrated 
N2H + 4-3 intensity (not shown) is found towards E33, while additional strong N2H + 4-3 emission 
is found towards E32 and the 850 /um continuum emission peak. A separate integrated intensity 
peak is found in northeast B2 (where we also find strong N2D + 3-2 emission). Overall, while 
the N2H + 4-3 emission is generally confined to areas where there is also continuum emission, the 
strongest N2H + 4-3 emission is found away from the brightest continuum contours. A comparison 
between Figure[3]and the N2H + 1-0 integrated intensity in Figuredb (greyscale) reveals a significant 
discrepancy between the locations of strong N2H + 1-0 and 4-3 emission in Oph B. We probe the 
possible causes of the offset, including the excitation and opacity of the lines, further below. 



4. Analysis 
4.1. Line fitting 

The rotational transitions of the linear molecular ions N2H + and N2D + contain hyperfine 
structure, which is dominated by the interactions between the m olecular electric field gradient 
and the electric quadrupole moments of the two nitrogen nuclei (jDore et al.l 12004! ) . Rotational 
transitions at higher J contain greater numbers of hyperfine components at smaller frequency 
intervals, i.e., the satellite components are found more closely spaced and can overlap significantly. 
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In particular, the N2H + 4-3 and the N2D + 3-2 tr ansitions contain 21 and 25 hyperfine c omponents, 



respectively, with relative line strengths > 0.001 (jPagani et al.1 12009 ; iGerin et al.ll2001 



For N2H + 4-3, we calculate the relative line strengths s j p 1 f^7' f^ p/ from the Einstein A 



coefficients Aj Fl f^,j> f[ f' an d frequencies vj Fl p^ji ^ p> given by Pagani et al. (j2009l ). using 



A, 



64vr 4 2 3 J 

^jl 1 V J F^F^J* F{F>Tp\ s J F X F-+J> F[. 



(2) 



where [i = 3.4 D is the permanent electric dipole moment, J is the upper rotational quantum number 
and [F] = IF + 1 accounts for the degeneracy of the hyperfine states. For the N2H" 1 " 4-3 transition, 
fourteen components at 8 distinct frequencies make up the main centroid emission, and lie within 
0.4 km s -1 of each other. At low optical depths the line profile does not differ substantially from 
a Gaussian profile, while the relative intensities of the satellite lines to the main line component 
remain small even at large opacities (~ 20% at r ~ 10). Nevertheless, the satellite components of 
the line are visible towards several pixels in Oph B (identified as locations A through C on Figure 
[3|), largely near the embedded protostar E33. In some locations the line structure is confused with 
a significant blue emission shoulder, which is not due to any hyp erfine structure, but w hich may 
be associated with an outflow stemming from either E33 or E32 (IKamazaki et al.ll2003l . White et 
al. 2010, in preparation). We were thus able to fit well the hyperfine structure towards only one 
location in the map, and constrain the excitation temperature, T ex , and opacity, r, towards two 
others. 

We fit the hyperfine structure with a custom Gaussian hyperfine structure fitting routine 
written in IDL, described in detail in Paper I, assuming equal excitation temperatures T ex for each 
hyperfine component. The line fitting routine does not take into account the significant overlap 
between hyperfine components in both the N2H + and N2D + transitions, but we find the total 
N2H + line opacities are still small where the satellite structure is seen (r ~ 2 — 3) and we therefore 
do not expect line overlap to introduce substantial errors. 

At low opacities, we cannot solve for r and T ex independently. Given a reasonable estimate 
for T ex , however, we may attempt to fit the line without visible satellite structure. In Paper II, we 
determined T ex for N2H + 1-0 emission across Oph B, and to first order can expect the same T ex to 
describe the excitation of all rotational levels of N2H + . This ass umption was found t o be accurate 
in a recent study of N2H + and N2D + emission in L1544, where ICaselli et al.1 (l2002bl ) were able to 
determine separately the excitation temperatures for N2H + 1-0, 3-2 and N2D + 2-1 and 3-2, and 
found they were consistent within uncertainties. We find, however, that while most of the data 
can be fit reasonably well using the N2H" 1 " 1-0 T ex values, visual inspection showed that a simple 
Gaussian fit the data equally well, or significantly better, in many cases, and the line opacities were 
not well-constrained. 



Because multiple components of similar line strength overlap to form the main N2H + 4-3 
component, a Gaussian fit will overestimate the line width, while the returned vlsr will be shifted 
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relative to a full HFS fit. The shift in vlsr, can be corrected (Avlsr = 0.036 km s _1 ), and at low 
opacities, the increase in Av is small (< 10 % for r < 1, up to ~ 25 % for r ~ 2). In the following 
sections, we thus discuss only the results of the Gaussian fitting of the N2H + 4-3 emission unless 
otherwise stated. 

No hyperfine structure was seen above the rms noise in N2D + 3-2 emission across Oph B. 
We are therefore unable to fit well the N2D^ 3-2 hyperfine structure solving for both t and T ex . 
Here, however, equating the N2D + T ex with that found for N2H + 1-0 emission resulted in better 
const raints on the opaci ty. The hyperfine fit, using the relative component strengths and velocities 



from iGerin et al.l (|200ll ). is thus more informative than a simple Gaussian fit, and we discuss the 



results from the full HFS fit in the following sections. 

To properly calculate the line width, Av, we subtracted in quadrature the resolution width, 



Av res (0.1 kms -1 for all species), from the observed line width, Av D b s , such that Av = J Av^ bs — Avj 
The uncertainties reported in the returned parameters are those determined by the fitting routine, 
and do not take the calibration uncertainty into account. The calibration uncertainty does not 
affect the uncertainties returned for v LSR or Av. The excitation temperature, column densities and 
fractional abundances discussed further below, however, are dependent on the amplitude of the line 
emission, and are thus affected by the absolute calibration uncertainty. 

The H2D + In — lio emission line has no hyperfine structure. We consequently fit the emission 
at each pixel with a single Gaussian (again in idl) to determine the line amplitude, vlsr an d Av 
across Oph B2. 



4.2. Centroid velocity and line widths 

Table [2] lists the mean, rms, minimum and maximum values of vlsr an d Av for N2H + 4- 
3, N 2 D+ 3-2, and H 2 D+ In - li , along with the same results from NH 3 (1,1) and N 2 H+ 1-0 
presented in Paper I and Paper II for comparison. 

We show in Figure [5] the N2D + v LSR and Av over Oph B2. We only plot values where the 
S/N ratio of the peak line intensity was > 5. The mean v LS r across B2 is 4.02 kms -1 with an rms 
variation of 0.19 km s -1 . Lower i> LSR values are found towards the western and northeastern regions 
of the Core, with the maximum v LS r = 4.44 kms -1 in the east towards the NH3 (1,1) clump B2-A7 
(identified in Paper I and labelled on Figure [5j coincident with the N2H + 1-0 clump B2-N6 from 
Paper II). The variation in v LSR does not suggest global rotation of B2. The N2D + line velocities 
shift between 3.7 kms -1 and 4.4 kms -1 over ~ 40", or ~ 4800 AU, to the west of E32. Nearby, the 
line f LSR shifts from 3.7 kms -1 to 4.2 kms -1 on either side of the western embedded protostar. 

The mean N2D + 3-2 line width Av = 0.53 km s _1 with an rms variation of 0.20 k ms -1 . In gen- 



eral, the widest line widths are found near E32 and the western protostar (Oph-emb5. lEnoch et al 



2009), to a maximum Av = 1.36 kms 1 in western B2 (but note the largest Av values are found 
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only ~ 45" to the northwest of some of the narrowest Av values). The locations of large line 
widths correspond well with the locations of dramatic vi sr shifts, described above. Away from the 
protostars, only small variations in the relatively narrow line widths are seen, notably towards the 
northeast and southwest, with Av ~ 0.3 km s -1 . Two objects are associ ated with localized narrow 



Av. An 850 fim continuum clump (162725-24273, iJprgensen et al.ll2008l ). shown on the image, was 
identified where the southwestern narrow line widths are found. In the east, B2-A7 NH3 is also 
associated with a Av minimum of 0.36 km s -1 . This is similar to the clump line width in N2H + 1-0 
emission at 18" resolution, but nearly a factor of two larger than found in N2H" 1 " 1-0 interferometer 
observations of the clump presented in Paper II. 

In Figured and Figure [6)3 we show the v LSR determined through Gaussian fitting of the H2D + 
In — lio an d N2H + 4-3 lines, respectively. Visual inspection of the v LSR values shows that the 
H2D + emission is similar kinematically to the N2D + 3-2 values in Figured, while a small shift in 
v LS r of ~ +0.2 km s -1 is seen in N2H + 4-3 emission towards Elias 32 and the 850 /im continuum 
emission peak and co-located 850 /um continuum emission peak (co-located with the continuum 
clump B2-MM8). The shift is small but significant given the otherwise excellent agreement in v LSR 
between the different tracers over B2 (see Table [2]). 

We show in Figure [6b and Figure EJl the Av of the H2D + and N2H + lines. The H^D" 1 " line 
widths are large, with a mean Av = 0.74 km s , and do not vary greatly over Oph B2 (Av rms 
= 0.19 kms^ 1 ). Slightly greater line widths are found in the western half of B2. Gaussian fits to 
the N2H + 4-3 data find the lines are slightly more narrow, with a mean Av = 0.69 km s -1 . As 
noted in §4.11 Gaussian fits likely overestimate the true N2H + line width by ~ 10 — 25 %. Wider 
N2H + lines are also found towards western Oph B2, while narrower lines are seen towards the 
submillimetre continuum emission peak and in B2's eastern half. N2H + lines of average width 
are found towards the embedded protostars, but become wider at pixels directly offset from the 
infrared sources. In particular, directly to the east of E32 and E33, the N2H + 4-3 line profiles 
show a significant blue shoulder with emission extending to v LSK ~ 2.5 km s -1 . Figure shows 
the result of a two-component Gaussian fit to the line and blue shoulder at location C in Figure El 



4.3. N2H + opacity and excitation temperature 

Comparison of the N2H + 1-0 and 4-3 line emission (see Figure and Figure [3]) shows that 
while in both cases emission is found over much of Oph B2, the integrated intensity distributions 
of the two transitions are quite different. It is thus unclear whether the excitation conditions for 
both transitions would be equal. 

Towards several pixels, hyperfine structure is visible. These locations are labeled A through C 
on Figure El and these spectra are again shown in Figure [H Using the full hyperfine component fit 
(overlaid on the spectrum), we were able to determine reasonably well the excitation temperature 
and opacity of the N2H + 4-3 line, T eXi 4_3 and T4-3 towards location A, with the resulting best fit 
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total line opacity T4_3 = 2.5 ± 1, with T ex ^_z = 8 ± 1 K. Both T eXj 4_3 and T4-3 are approximately 
consistent with values found in Paper II for N2H + 1-0 emission at the same location, where 
T ex ,i~o = 7.3 ± 0.1K and ti_o = 3.1 ± 0.6. The brightest N2H + 4-3 line emission, however, 
and additionally the maximum integrated intensity, is located towards E33 (location B in Figure 
[3l directly south of A). The intensity ratio of the main N2H + 4-3 component with an apparent 
red satellite (line strength only 2x the data rms noise level, but at the correct velocity and with 
a width that matches the main component, see Figure HJj) suggests that the N2H + 4-3 opacity 
towards E33 remains moderately optically thick, T4_3 ~ 3, and the corresponding T^j 4—3 ~ 8. The 
expected line profile given these values is overlaid on the spectrum in Figure [U The integrated 
N2H + 1-0 emission decreases significantly towards E33 (see Figure [lb and Paper II). Additionally, 
the emission is confused by multiple velocity components which preclude robust fits to the 1-0 
spectra, so we do not have a good estimate of the N2H + 1-0 T ex at this location. At location D 
in Figure [3l satellite emission is again seen, but here the hyperfine fits are confused by significant 
blue shoulder emission. The best-fit results give T ex <~ 7.5 K — 10 K, with t ~ 1 — 3. We overlay the 
expected line profile for T ex = 8.5 K and r = 2 on spectrum D in Figure HI Elsewhere, no hyperfine 
structure is visible in the N2H + 4-3 line. 

For opacity ratios between the N2H + 4-3 and N2H + 1-0 lines, t^-^/ti-o ~ 0.3 and similar line 
widths, the derived N2H + column density, A r (N2H + ), calculated using either transition following 



Caselli et al.l (|2002d ) (Equation Al), is equal if both T ea;j i_o and T ex ^_z ~ 7K. At lower opacity 



ratios, the required T ex also decreases. On average, the N2H + 1-0 and 4-3 line widths agree, 
however we note that since we cannot fit well the hyperfine structure, the N2H + 4-3 widths are 
an overestimate of the true value by ~ 10 — 25 %. This uncertainty propagates linearly into the 
column density calculation. With this caveat, an excitation temperature T ex ~ 7K agrees well with 
the mean T ex determined from HFS fitting of the JN^H^ 1-0 line in Paper II. We have no column 
density estimate towards E33 from N2H + 1-0 emission. Near E32 and E33, however, the greater 
relative N2H + 4-3 opacity requires either a higher common T ex ~ 10 — 11 K to obtain similar 
A^(N2H + ) estimates, which does not agree with our N2H + 1-0 T ex measurement, or the excitation 
of the N2H + 4-3 line must be greater than that of the 1-0 line by ~ 2 — 3 K. This moderate increase 
in excitation temperature of the higher excitation N2H + line is not ruled out by our data. Detailed 
modeling of the N2H + emission near the embedded protostars using radiative transfer codes would 
help disentangle the excitation and opacity effects. 



4.4. Column density and fractional abundance 



We next discuss the column densities of N2D + and H^D" 1 " derived from the observed line 
emission. The N2H + column density was calculated in Paper II. 
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44-1- N 2 D+ 

We calculate the to tal N2D" 1 " column de nsity, A^ t(N2D + ), from the integrated intensity of the 



3-2 transition following ICaselli et al.l (|2002al ) for optically thin emission: 



Nt t = 87tW 9u 1 1 Qrot(T ex ) 

\ 3 A gi J u (T ex ) - J u (T bg ) 1 - exp(-hv/kT ex ) giexp(-Ei/kT ex ) 

where W is the line integrated intensity, A is the Einstein spontaneous emission coefficient, g u 
and gi are the statistical weights of the upper and lower states, respectively, and the equivalent 
Rayleigh-Jeans excitation and background temperatures are given by J u (T ex ) and J u (Tb g ). The 
rotational energy Ej = J{J + l)hB for linear molecules, where the rotation quan tum number 



J = 2 for the lower energy state and B = 38554.717 MHz is the rotational constant (pastry et al 



1981 



). The partition function Q ro t{T ex ) = ]CjLo(2^+l) ex P( — Ej/kT ex ), and we use the integrated 
value, Qrot(T ex ) = kT ex / (hB) + 1/3. We assume that the excitation temperature T ex is the same 
for all rotational levels, and use the fitted N2H + 1-0 T ex from Paper II as an estimate of the 
excitation temperature of the N2D + line at each pixel (as discussed in §4.ip . All other parameters 
remaining equal, at the average T ex of the N2H^ 1-0 line, a decrease in T ex by 2 K would increase 
the resulting N2D + column density by a factor < 2, while an increase in T ex by 2K would decrease 
the resulting N2D + column density by < 0.5. We show iV(N2D + ) across Oph B2 in Figure [7^,. 

We then calculate the fractional N2D + abundance per pixel, 



A(N 2 D+) = iV(N 2 D+)/iV(Il2) (4) 

where the H2 column density, N(H.2), is calculated from 850 /um continuum emission. The 15" 
data were convolved to a final 18" FWHM resolution to match the convolved N2D + data. The 
H2 column density is then given by iV(H2) = S u /[Q m fim H K u B u (T c i)], where S u is the 850 ^m flux 
density, Q m is the main-beam solid angle, \i = 2.33 is the mean molecular weight, m H is the mass of 
hydrogen, k v is the dust opacity per unit mass at 850 fim, where we take k v = 0.018 cm 2 g _1 , and 
B u {Td) is the Planck function at the dust temperature, Tj. As in Paper I and Paper II, we expect 
the gas and dust to be well-coupled in Oph B2, and set = Tk, the gas kinetic temperature 
determined from NH3 (1,1) and (2,2) inversion line ratios. We refer the reader to the detailed 
discussion of the uncertainties in the iV(H2) values presented in Paper II, but note that the derived 
H2 column densities have uncertainties of factors of a few. Additionally, due to the negative artifacts 
introduced into the data from the observational chopping technique (again discussed in detail in 
Paper II), we limit the analysis to pixels where S u > 0.1 Jy beam -1 , though the rms noise level 
of the continuum map is ~ 0.03 Jy beam -1 . For a dust temperature = 15 K, this flux level 
corresponds to iV(H2) ~ 6 x 10 21 cm -2 . 

The column density distribution is not similar to the the integrated intensity map shown in 
Figure [Th. Figured shows several small peaks of A^(N2D + ) in B2 which largely correspond with 
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local minima in the excitation temperature (T ex ~ 5 K compared with the mean T ex ~ 7 K , where 
the uncertainty in T ex is ~ 0.4 — 0.6 K). The 7V(N2D + ) peaks remain within the continuum contours 
but avoid maxima of continuum emission. We find a mean iV(N2D + ) = 1.8 x 10 11 cm -2 in Oph 
B2 with an rms variation of 1.3 x 10 11 cm -2 . The maximum 7V(N2D + ) = 7.1 x 10 11 cm~ 2 is found 
towards the northeast. If we calculate A r (N2D + ) assuming a constant T ex = 7K, the distribution 
follows the integrated intensity but the mean value and variation remain the same. 

The fractional N2D + abundance (not shown) follows a similar distribution to the N2D + column 
density. We find a mean X(N2D + ) = 8.2 x 10 -12 with an rms variation of 5.1 x 10 -12 . The maximum 
X(N2D + ) = 2.9 x 10~ n is again found towards northeast Oph B2. Within the continuum emission 
contours, the lowest measured N2D + abundances (X(N2D + ) ~ 2 x 10~ 12 ) are found towards the 
continuum emission peak and B2-MM8 clump (emission towards the Elias 33 protostar was not 
above our S/N limit for analysis). 



44.2. H 2 D + 

To calculate the H2D + column density, iV(H2D + ), we first estimate the line opacity r from 
the H2D + In — lio line emission using 



In 



Tmb 



J {Tex) — J{Tbg 



(5) 



where we must assume a priori an excitation temperature T ex for the transition. An upper limit 
is given by the gas kinetic temperature Tk- For the pixels where the H2D + S/N ratio > 5, we 
found in Paper I a mean Tk = 14.5 K with an rms variation of 1.1 K (slightly less than the mean 
Tk = 15 K found for the entire Oph B Core in Paper I), with a minimum Tk < 13 K towards 
the highest 7V(H2) values. We note, however, that the critical density at which NH3 is excited 
(n cr *~ 10^ ^ cm ^) is lower than that of the H^D^ transition by an order of magnitude, and the 
NH3 temperatures, or the variation in Tk, may not accurately reflect the conditions where H2D + 
is excited. 



In their recent H2D + study of several starless and protostellar cores, ICaselli et al.l (|2008l ) found 



that the H2D + T ex was similar to or slightly less than the core Tjc values, which were estimated 
through a variety of methods (e.g., continuum observations, NH3 inversion observations as for B2, 
and multiple transition observations of molecules such as HCO + ). Six of the seven protostellar cores 
in their sample had excitation temperatures which ranged from 9K to 14 K and similar densities 
to th ose calculated for B2 in Paper II. Based on these values, and to better compare our results 
Caselli et all we use a constant excitation temperature T ex = 12 K in Equation [5j Given 
: 12 K, we fin d a mean r = 0.13, with a minimum r = 0.05 and maximum r = 0.28, which 



with 
T - 

- 1 ex 

agree well with the ICaselli et al.1 results in protostellar cores. 
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The total column density of ortho-H^D" 1 " is then given bv lVastel et al.l ((2006) 



iV(ortho-H 2 D H 



8tt Qrot(T e x) exp(E u /kT e . 



A 3 A 



ul 



9u 



exp(hv/kT e , 



rdv 



(6) 



where g u = 9, E u /k = 17.8 K, and A u i = 1.08 x 10 4 s 1 for the li 1 — lm transition (iRamanlal Tennyson 
2004 ). The integral f rdv = \^/tt / {\n.(2)T&.v. W e calculate Q rot by reducing the H 2 D + level 



structure to a 2-level system following ICaselli et al.l Since the energy of the first excited state above 
ground is E/k = 17.8 K and that of the second excited state is E/k = 110 K, we expect this approx- 
imation to be valid in Oph B2 given the low temperatures determined in Paper I. The partition 
function depends on T ex , which we have estimated to be 12 K. In Figure [HJ we plot the variation 
in r, Q ro t and iV(ortho-H2D + ) with T ex given a line Tmb = 0.5 K and Av = 0.7 km s^ 1 (typical 
values in B2), and have listed the returned values for T ex between 7K and 15 K in Table [3l While 
the variation in returned parameters is fairly large with small T ex changes at low T ex values (i.e., 
factors ~ 2 or more around T ex ~ 7K), at higher T ex the values vary less. For example, an increased 
T ex = 14 K rather than T ex = 12 K would result in a decreased iV(ortho-H2D + ) by ~ 20%, while a 
decreased T ex = 10 K would result in an increased 7V(ortho-H2D + ) by ~ 30%. We thus assume an 
uncertainty of ~ 25% in the iV(H2D + ) values reported here. 

To calculate the total H2D + column density in B2, an estimate of the ortho- to para-H2D + 
(o/p-H_2D + ) ratio is needed. The ground state para-H^D" 1 " transition occurs at a frequency of 
1.37 THz; this line is difficult to observe from the ground and may not be excited at the low 
temperatures found in star-forming cores. Without a direct measurement of the population of the 
para state, we look to chemical models to pred ict the o/p-H 2 p + rati o, which show that the o/p- 
H2D + ratio is directly dependent on o/p-H2. IWalmslev et al.l (|2004l ) find that the o/p-H2 ratio 



increases to a steady-state value of ~ 5 x 10 -4 at T = 10 K and n = 10 6 cm -3 , similar to the 
density and temperature determined for Oph B2 in Paper I and Pa per II. The same m odels predict 
o/p-H2D + ~ 0.2 — 0.3. Using Bonnor-Ebert sphere core models, ISipila et al.l (|2010l ) also predict 
o/p-H2D + ratios < 0.5 at T ~ 10 — 15 K in steady-state. Their models show that the o/p-H 2 D + 
ratio can dramatically increase to values ~ 2 — 3 at T < 10 K, but these temperatures are unlikely 
based on our analysis of Oph B2. 

It is not clear, however, whether steady-state models accurately represent conditions in star 
forming cores like B2. The dynamical timescale of a collapsing core can be much shorter than some 
of the chemic al timescales, leadin g to large depar tures of free-f all model predictions from steady- 
state models (jFlower et al.ll2006l ). For example, iFlower et al.l show that the steady-state model 
underpredicts the o/p ratio by up to an order of magnitude compared with a free-fall collapse model 
at gas densities n < 10 6 cm -3 (the models are more consistent at higher densities). Consequently, 
it may not be unreasonable to expect variations in the o/p-H_2D + ratio across B2 as it evolves. 
Given this uncertainty in the o/p-H 2 D + ratio, we limit discussion to iV(ortho-H2D + ), but note 
that the total iV(H 2 D + ) is likely greater than A r (ortho-H 2 D + ) by a factor of a few. 



-15- 



We show in Figure [9] the resulting iV(ortho-H 2 D + ) distribution across Oph B2. Including 
only pixels where the S/N > 5, we find the mean 7V(ortho-H 2 D + ) = 1.4 x 10 13 cm -2 with an 
rms variation of 0.7 x 10 13 cm -2 over Oph B2. To this sensitivity limit, we find a minimum 
iV(ortho-H 2 D+) = 4.1 x 10 12 cm~ 2 . The maximum iV(ortho-H 2 D + ) = 3.3 x 10 13 cm" 2 is found 
~ 30" to the northwest of the 850 /um continuum peak. The largest ortho-H 2 D + column densities 
are found in the western half of Oph B2, with moderate iV(ortho-H 2 D + ) values extending to the 
northeast. 

We next calculate X(ortho-H 2 D + ) as described for N 2 D + above, and show in Figure [9)3 the 
A(ortho-H 2 D+) distribut ion across B2. Similarly to the column density distribution, the largest 
H 2 D + abundances, A(ortho-H 2 D + ) ~ 4 — 5 x 10~ 10 , are found in the west. A local X(ortho-H 2 D + ) 
minimum of ~ 2 x 10~ 10 is found towards the central continuum emission peak, with similar 
abundances towards the northeast. 



4.5. Using H 2 D+, N 2 H+ and N 2 D+ to determine T K 

The molecular weight of H 2 D + (m mo i = 4.01 m# ) is much less than that of N 2 H + or N 2 D + 
(29.01 m# and 30.02 m^, respectively). If we assume that the H 2 D + In — lio, N 2 H + 4-3 or 
N 2 D + 3-2 transitions are excited in the same environment within Oph B2, we would expect each 
to trace the same core motions, and be described by the same kinetic temperature Tk- This co- 
location would imply that the non-thermal line widths for the two species should be equal, where 
°ivT = a obs ~ a th ■ Given equal Tk, the H 2 D + thermal line width, a t h = (fcTfe/wimoj) > wu l 
be measurably larger than that of N 2 H + due to its smaller molecular weight. At Tk = 12 K, for 
example, ath = 0.06kms _1 for N 2 D + and 0.16kms _1 for H 2 D + . We can use this difference in 
thermal line widths to determine the kinetic gas temperature by equating the non-thermal line 
widths and solving for Tk (e.g., in the case of N 2 H + ): 



T K 



1 



1 



m. 



[a 2 6s (H 2 D+)-a 2 )S (N 2 H+)] 



(7) 



Mean temperatures of 9.2 ± 0.2 K, in good agreement with tem peratures derived from NH 



alone, were derived in this manner towards a number of dense cores bv lFuller &: Mverd (|1993l ) using 
the molecular species HC3N and NH3. The largest source of error in the derived Tk likely arises 
from the assumptions listed above rather than from the propagation of the small uncertainties on 
the observed line widths. 

Overall, we expect the comparison of H 2 D + and N 2 D + to be most accurate, given the similar 
critical densities of the transitions, the fact that N 2 D + must form from H 2 D + , and the similar 
distribution of integrated intensity (in contrast with the poor correlation between H 2 D + and N 2 H + 
4-3 emission; see Figures QJ, and [2]). The different distributions of H 2 D + and N 2 H + 4-3 in Oph B2 
suggest that these two species are not excited in the same environment in the core. Additionally, 
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since our estimates of the N2H + line widths are not accurate (given the unresolved hyperfine 
structure), we will focus only on N2D + for this analysis. 

We show in Figure [101 the temperatures calculated based on Equation [7J for H 2 D + and N 2 D + 
3-2. Only pixels where the signal-to-noise ratio of both lines was > 5 are shown. The derived 
Tr values span a large range, from unphysical negative values (~ — 11K) to unrealistically large 
values (~ 75 K). Using the mean Av values of H2D + and N 2 D + from Table [2j we calculate a mean 
Tk = 16 K, similar to that determined through NH3 analysis. Multiple pixels with Tk ~ 12 K 
are found surrounding the central continuum emission peak. Negative temperatures, however, are 
found associated with the continuum peak and associated clump B2-MM8. 

From Equation we calculate negative Tk values when the observed H2D + line widths are 
less than those of N2D + 3-2. We show in Figure [TT] (bottom) Av derived for H2D + and N 2 D + 
emission (data values in 15" pixels plotted), along with the spectra and relative fits towards B2- 
MM8 (top). Over most of B2, N 2 D+ line widths are narrower than H2D + , and thus Tk > OK. The 
negative temperatures derive from two adjacent pixels towards B2-MM8 where the N2D + Av > Av 
of H2D + . Figure [TT1 shows non-Gaussianity in the spectra of both species towards B2-MM8, which 
complicates the calculation of a single line width to describe the emission. 

Why would we find greater N2D + Av than H2D + Av only towards the central continuum 
peak? One possibility is suggested by the spectra in Figure \TT\ where additional emission in the 
H2D + spectrum is clearly seen blue-shifted relative to the main Gaussian peak. For the H2D + 
line, the best-fit Gaussian does not include the blue line shoulder. Slight blue asymmetry is also 
visible in the N2D + line, but the distinction between the main component and the shoulder is less 
obvious. It is possible that the shoulder emission of N2D + 3-2 is spread over a larger range in v LSR 
due to the contribution of multiple hyperfine components to the observed emission (which would 
each contain emission at the central line velocity and at the shoulder velocity), broadening the line 
overall. The best fit to the N2D + line thus cannot distinguish well between the main and shoulder 
emission, resulting in an artificially broad Av. 



5. Discussion 

5.1. Line widths and density 

The mean line widths listed in Table [2] show that on average, the observed Av for a given 
molecular line decreases with increasing critical line density, suggesting that the gas in B2 becomes 
more quiescent at higher densities (omitting ^H -1- 4-3, where Av is an overestimate). For an 
average gas temperature Tk = 15 K, the resulting non-thermal line contribution decreases from 
<t nt = 0.36 km s -1 for NH3 (1,1) emission to <r NT = 0.22 km s -1 for N2D + 3-2 emission (for 
n cr ~ 10 3-4 cm -3 to n cr ~ 8 x 10 5 cm -3 , respectively). Comparing with the thermal sound speed 
at 15 K, where c s = y / kTk / (fJ-m H ) , we find an average <7 NT / c s = 1.6 for NH3 and a NT / c s = 1.1 
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for N2D + 3-2 emission. This comparison shows that despite the decrease in non-thermal motions, 
the non-thermal motions are equal in magnitude with the sound speed, even at high density. A 
lower gas temperature would result in even larger magnitudes of <t nt , while significantly higher 
temperatures are unlikely. 

T his decrease in non-thermal motions at higher densi ties has been observed in starless cores 



(e.g., iPon et alj|2009l : iTafalla et al.ll2004j : lLada et al.ll2003l ). Frequently in starless cores, however 



the narrowest observed line widths are consistent with being due to nearly pure thermal motions, 
in contrast to the results found here. As a NT values decrease, the ability for turbulent pressure 
to support the core against collapse is reduced, which impacts the stability of the core against its 
self-gravity. The narrow <r NT widths found in starle ss cores are often expl ained as the result of the 
dissipation of turbulent motions in the core centre (IGoodman et al.lll998l ). If the non-thermal line 



widths in B2 are due to turbulent motions, the corresponding turbulent pressure (Pnt 
is equal to or greater than the thermal pressure (Pp = nkT). 



mnaz 



In Paper I, we discussed the source of wide NH3 line widths in Oph B, and concluded that if 
the relatively large <r NT values are caused by turbulent motions, they could not be primordial (i.e., 
inherited from the parent cloud). This is because the dampin g timescale of such motion s (with no 
external driving force) is approximately the dynamical time (IMac Low &: Klessenl 12004 ). and the 
existence of protostars associated with B2 implies that B2 has likely existed, with high density gas, 
for at least a dynamical time. 

The wide lines in Oph B2 may, however, be due to more organized motions such as infall 
or outflow. None of the observed molecular lines are optically thick, and we therefore do not 
see the self-absorbed, asymmetric line profiles frequently used to infer the infall of gas within a 
core. Infall si gnatures were obser ved in CS, H2CO or HCO + emissi on towards several lo cations in 
eastern B2 by I Andre et al.l (120071 ). but not in central or western B2. iGurney et al.l (|2008l ) observed 
complicated CO (and CO isotopologues) line structures towards the B2 continuum peak, but were 
unable to distinguish infall conclusively. The analysis of line profiles is complicated by a substantial 
outflow discovered in CO emission and emanating from one of the two embedded YSOs in B2 
( Kamazaki et al. 2003 ) which extends over 10' (~ 0.4 pc) in siz^E Upcoming results of CO, 13 CO 
and C 18 3-2 observations i n Ophiuchus (White et al. 20 10, in preparation) from the JCMT Gould 
Belt Legacy Survey (GBLS, IWard-Thompson et al.ll2007l ) will be better able to study the infall and 
outflow motions in Oph B. In conjunction with these data, the process of how the outflow motions 
affect the dense gas may be studied. There is clear evidence, however, for protostellar influence on 
the N2D + line widths shown in Figure [5b, as discussed in §4.21 Additionally, the blue shoulder 
emission seen in N2H + and N2D + east of E32 and E33 coincides with the eastern CO lobe of the 
outflow, and is further evidence for the significant impact of the outflow on the dense gas. 



We also noted in §4.21 that there is an offset in v LSR of ~ 0.2 km s 1 between N2H + 4-3 and 
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H2D + emission towards B2-MM8 and the nearby embedded protostar Elias 32, with the N2H + 
emission redshifted relative to the H2D + emission. Over the rest of Oph B2, the correlation 
between the v hSR of the two lines is good, suggesting the offset is real and localized. Comparison 
with N2D" 1 " emission at the same location shows that the v LSR of the H2D + and N2D + agree well, 
and it is consequently the N2H + emission which is offset in velocity from the rest of the core gas. 
N2H + 4-3 has a critical density approximately an order of magnitude greater than the H^D" 1 " and 
N2D + lines discussed here. The offset in w LSR may thus be probing a shift in the gas motions at 
extremely high densities in the core interior. 



5.2. Trends in the deuterium fractionation 



5.2.1. The N 2 D + /N 2 H + ratio 



The deuterium fractionation, Rd, can be defined as the ratio between the column den- 
sities of a deuterated molecule and its hydrogen-bearing counterpart. Here, we define Rd = 
iV(N 2 D+)/iV(N2H + ), and calculate R D in Oph B2 from iV(N 2 D + ) calculated in gXTjand iV(N 2 H + ) 
calculated in Paper II. Figure [7b shows Rd across B2. The distribution is similar to the N2D" 1 " col- 
umn density in Figure [Th. The largest values of Rd, to a maximum Rd = 0.16, are found towards 
northeast B2, while towards the smaller scale maxima, we find Rd ~ 0.1 or less. We estimate 
we are sensitive to Rd ~ 0.01 — 0.015 based on the rms noise in the N2D + integrated intensity 
map (3cr ~ 0.2Kkms _1 in the off-B2 pixels), and the mean N2H + column density over Oph B2, 
(A r (N2H + )) = 5.9 x 10 12 cm -2 , calculated in Paper II. We find a moderate mean Rd = 0.03 with an 
rms variation of 0.01 towards B2, comparabl e to results in Ori B9 t owards t wo N2H" 1 " emi s sion p eaks 
and a protostellar source (R D = 0.03-0.04, llVIiettinen et al.l Eo09^. While IPagani et al.l (l2009h find 



substantially greater central Rd values through modelling of the starless core L183 (Rd ~ 0.7 at 
the core centre) which decrease to ~ 0.06 at ~ 5000 AU, it is not clear how to compar e their model 
of Rd as a function of core radius to our observed column density ratio. iFontani et al.1 (|2006l ) found 
an average Rd ~ 0.01 towards the high mass star forming region IRAS 05 345+3157, which at high 
resolution resolved into t wo N2D" 1 " co n densa tions each with Rd = 0.11 (IFontani et alJ 120081 ). In 
a study of starless cores, ICrapsi et al.1 (120051 ) find a range in Rd between a low er limit of 0.02 to 
a maximum of 0.44 (which was found towards the Oph D Core). ICrapsi et al.1 also note that Rd 



values > 0.1 are generally only found for cores with iV(N2H + ) > 10 13 cm -2 , which is approximately 
the maxi mum A f (N2H + ) in Oph B2. A similar range in Rd was found towards protostellar, Class 
sources (jEmprechtinger et al.l 12009). On average, then, our results agree with previous results. 
Given our map of Oph B2, however, we can also probe the variation of Rd across the core. 

In Figure [T2l we show Rd as a function of increasing iV(H2) (Figure fT2b). and also as a function 
of increasing distance (in projection) from the nearest embedded protostar (Figure Eb). Each data 
point represents an 18" pixel (to match the beam FWHM) and pixels with S/N < 5 in either N2H + 
or N2D + emission are omitted. No clear trend is present in the deuterium fractionation with iV(H2). 
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We find, however, that Rd increases at greater projected distances from the embedded sources, as 
might be expected from the N2D + integrated intensity distribution. We fit a linear trend to the 
data, where Rd = (—0.007 ± 0.002) + (0.37 ± 0.06)d[pc], but note that there is large scatter in 
the relation. An analysis of the individual N2H + and N2D + abundances shows that the increasing 
Rd trend with projected protostellar distance is dominated by a similar increase in A(N2D + ) with 
projected protostellar distance, while A(N2H + ) is approximately constant over the range plotted, 
shown in Figure \l2ti. 



5.2.2. H 2 D+ 

We next look at the distribution of the ortho-H^D" 1 " abundance over Oph B2. In Figure [T2l 
we show X(ortho-H2D + ) as a function of iV(H2) (Figure fT2e) , and also as a function of increasing 
projected protostellar distance (Figure [T2T). We show data in 15" pixels, and omit pixels with 
S/N < 5. We find no trend in the fractional ortho-H^D" 1 " abundance with H2 column density (an 
increasing trend is found, but not shown, in the ortho-H2D + column density, iV(ortho-H2D + ), 
with iV(H2)). The X(ortho-H2D + ) vs. iV(H2) relationship is not entirely a scatter plot, as 
the ortho-H2D + abundance tends to increase with projected protostellar distance to a maximum 
X(ortho-H 2 D + ) = 4.9 x 10" 10 at d = 0.035 pc, with A(ortho-H 2 D+) ~ 1.8 x 10~ 10 at larger dis- 
tances. All data points at d > 0.035 pc correspond to the emission found in northeast B2, while the 
rest show emission closer to the central continuum emission peak. If we fit a linear relationship to the 
data, omitting the northeast B2 points, we find A(ortho-H 2 D+) x 10~ 10 = (1.2±0.5)+(79±23)d [pc] . 

There is also no clear trend in iV(ortho-H2D + ) with Rd where both H^D" 1 " and N2D + are 
strongly detected across Oph B2, although both species are found with relatively large column 
densities north and west of the continuum peak. Given the small-scale structure found in the 
iV(N2D + ) distribution in B2, the spatial resolution and sampling of our H?D + observ ations may 



not be high enough to discern a correlation. We note, however, that lCaselli et al.l (|2008l ) also found 
no clear correlation between Rd and iV(ortho-H2D + ) in their sample of starless and protostellar 
cores (with the caveat that some Rd values were derived using column density ratios of deuterated 
NH3 and deuterated H2CO with their undeuterated counterparts rather than iV(N2D + )/iV(N2H + )). 
In particular, high values of Rd and relatively low A r (ortho-H2D + ) were found towards two objects 
in Ophiuchus (Oph D and 16293E). 



5.3. What is affecting the deuterium fractionation in Oph B2? 

The largest deuterium fractionation values in Oph B2 are not found towards the continuum 
emission peak, but instead avoid the continuum peak and nearby protostars. The results of the 
previous section indicate that the distance to the nearest protostar is the dominant parameter im- 
pacting the deuterium fractionation in Oph B2. We next discuss several mechanisms by which the 
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embedded sources are most likely to impact the deuterium fractionation distribution in B2, includ- 
ing temperature variation in the gas, liberation of CO from dust grains, the impact of protostellar 
outflows, and possible increase in the ionization fraction due to x-rays from the protostars. 



5.3.1. Temperature and CO depletion 

Temperature can have a substantial impact on the deuterium fractionation in a dense core. As 
stated in £ )4.4,H we expect the gas and dust to be well-coupled at the densities of Oph B2 such that 
Td = T gas , but specify in the following which of the gas or dust temperature is important in the 
following mechanisms. First, at T gas > 20 K, the H2D + reaction in Equation [1] can proceed both 
forwards and backwards, resulting in no net increase in H2D + , and consequently other deuterated 
species formed via reactions with H2D + . Second, the CO which was deposited (adsorbed) as an icy 
mantle onto dust grains while at low temperatures will evaporate (desorb) back into the gas phase 
if the dust grains are heated, providing a destruction mechanism for H2D + and again interrupting 
the deuteration chain. 

The CO depletion factor, /d(CO), is defined as the canonical CO abundance, X can (CO), 
divided by the CO abundance, X{CO) = iV(CO)/7V(H2), in the observed region: 



/d(CO) 



X C an (CO) 

X(CO) 



(8) 



such that fr> values are low where observed CO abundances are high. Low CO depletion factor s 
can dramatically lower the H2D + abundance at gas temperatures T < 15 K (jCaselli et al.ll2008l ). 
Currently no estimates of fjj exist in B2 due to a lack of data at spatial resolutions matching 
the available continuum data. The CO observations of the JCMT GBLS should enable the first 
calculations of fr> at complementary spatial resolution, and allow further analysis of the relationship 
between CO and the deuterium distribution in B2. The extended distribution of B\2D + and N2D + 
across Oph B2 implies, however, that extensive depletion of CO has occurred. At densities of 
n = 10 6 cm -3 , CO adsorption and desorption occur equally frequently when Td = 18 K ( Visser et al 



2009). At Td > 18 K, CO will thus again become significantly more abundant in the gas phase, 



while at T gas > 20 K, enhanced deuterium fractionation via Equation Q] will cease to operate. 

An obvious potential source of heating is the presence of embedded protostars within Oph B2. 
In a ID dynamic /c hemical model of an isolated dense core which collapses to form a protostar, 
Aikawa et al. (|2008l . see their Figure 2) find temperatures greater than 18 K can be found at radii 
~ 5000 AU (~ 0.02 pc) from the protostar at ~ 10 5 yr after its formation. As a function of total 
luminosity of the core, L core , however, their model predicts a CO sublimation radius of < 1000 AU 
(0.005 pc) for L 

core ~ 0.5 — 1 Lq. In a more specific example, IStamatellos et al.l (|2007l ) model the 
heating effect of an embedded Class I protostar (Lj, i = 10 Lq) on the surrounding gas in the Oph 
A Core, and find only a small temperature increase (< 5K), limited to the immediate vicinity of 
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the protostar, in agreement with lAikawa et al.l ((2008) 



Bot h YSOs in central B2 have been classified as embedded (jj0rgensen et al.ll2008l ). Class I 
sources (jEnoch et al.l 120071 ). The estimated bolometric luminosities, Lh„ i ~ 0.5 L and ~ 1 L 



for E32 and E33, respectivel y ([Enoch et al 



2007; 



van Kempen et al.l l2009). are much smaller than 



used in the heating model by Stamatellos et al. ( 20071 ) . The bolometric temperatures, ~ 300 K 
and ~ 500 K, suggest both are in the later stages of the Class I phase, as do their negative infrared 
spectral indices {am = —0.03 and —0.12). Sources with —0.3 < aiR < are sometimes classified 
as 'flat-spectrum' sou rces, and are thought to be an intermediate stage between Class I and Class 
II. ICrapsi et al.l (|2008l ) suggest that most or all flat-spectrum sources could be Class II objects seen 



edge-on. Based on these results, E32 and E33 are not likely to be deeply embedded within Oph B, 
and have ages ~ 5 x 10 5 yr or greater based on lifetime estimates of the different protostellar stages 



(|Evans et all 120091 ) . 



The low luminosities of E32 and E33 suggest that they should not have a large impact on the 
temperature of B2, in agreement with our temperature results in Paper I which showed very little 
variation in the gas temperature over Oph B. Figure [T2l shows that we find no trend in the kinetic 
temperature Tk, determined using NH3 inversion transitions in Paper I, with protostellar distance 
above the ~ 1 — 2K uncertainties in the data. The derived temperatures are near the required 
values, however, and the NH3 temperatures may not be indicative of the conditions of the denser 
gas where we expect most of the deuterated species to exist. Depletion of H2D + and N2D + via CO 
sublimation would thus require an increasing temperature gradient towards higher densities near 
the protostars, which is not unreasonable. 

Since CO also destroys N2H + , significant evaporation of CO from dust grains should also be 
observable in the N2H + distribution. Accordingly, in Figure [lb we see a decrease in the integrated 
N2H + 1-0 intensity towards E33. No good fit to the N2H + spectra was found at this location, 
however, due to complicated line structures, such that we are not able to determine the N2H + 
column density or abundance. Figure fT2l i shows no trend on larger scales of X(N2H + ) with 
protostellar proximity. The N2H + 4-3 distribution shows N2H + remains in the gas towards E33, 
but an analysis of the excitation of the line would be needed to determine the N2H + column density 
from these data. 

In summary, heating by an embedded protostar can produce slightly greater gas temperatures 
and evaporate CO from the surface of dust grains, which would decrease the deuterium fractionation 
in the gas. Since we do not see evidence of a decrease in X(N2H + ) towards the protostars, the 
evaporation of CO from dust grains is likely not the major mechanism for the decrease in abundance 
of deuterated species. While we also don't see evidence of increased gas temperatures near the 
protostars in the NH3 data, these temperatures reflect average values along the line-of-sight, and 
do not rule out increased temperatures at higher densities near the protostars. Given the strong 
dependence of increased deuterium fractionation on temperature, we expect this mechanism is likely 
very important in B2. 
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5.3.2. The Oph B2 outflow 

The deuterium fractionation may be impacted by the known, extensive CO outflow in B2 
described in £ 15.11 It appears that the dense gas in eastern Oph B2 is affected by the outflow, given 
the larger line widths seen towards the protostars in high density tracers, changes in vt,sr, and the 



blue shoulder emission which is coincident with the eastern lobe of the outflow (Kamazaki et al 



20031 ) (although large line widths could also indicate the presence of infall motions). We do not find 
evidence of outflow interaction with the dense gas in northwest Oph B2, suggesting that the outflow 
may be in front of or behind the bulk of the core gas. Very localized heating of the surrounding gas 
by the protostar could evaporate CO from the dust grains. The outflow, then containing entrained 
CO, could mix the CO with the dense core gas, destroying the deuterated species without an 
increase in the gas temperature on larger scales. Alternatively, shocks associated with the outflow 
could raise gas temperatures and reintroduce enough CO to the gas phase to affect the deuterium 
fractionation. Further analysis of the outflow, including observations of shock tracers, is needed to 
determine whether this scenario is reasonable, and can examine the impact of low- mass protostellar 
outflows on dense gas chemistry. As discussed above, however, we do not find evidence for variation 
in the N2H + abundance as would be expected if destruction due to gas-phase CO is the dominant 
mechanism for decreasing the abundance of deuterated species relative to their non-deuterated 
counterparts. 



5.3.3. Ionization fraction and electron abundance x e 



Variations in the ionization fraction in the dense gas can also affect the deuterium fractionation, 
in particular in regions of high CO depletion. The electron abundance, x e = n e /n K2 , is the ratio 



of the electron number density 



with the number density of molecular hydrogen, 



n. 



and 



ae reduced b\ 


r an increase 


(Caselli et al. 


2008). The 



in x e , which in turn increases the rate of dissociative recombination (ICaselli et al. 
abundances of deuterated molecules can thus be used as a probe of x e in dense cores. At high 
extinction towards starles s cores, the c reation of ions is expected to be dominated by the cosmic 



— 1/2 

ray flux, with x e oc n H2 (|McKedll989l ) 



Studies investi gating x e in dense core s using the DCO + /HCO + ratio have found average line- 
of-sight x e ~ 10 -7 (jBergin &: Tafallall2007l . and references therein). In regions of high depletion, the 
column density ratio of N2D + with N2H + is expected to be a better probe of x e than the ratio of 
DCO + with HCO + . Since is expected to be the major molecular ion in cores strongly depleted 
in heavy elements, we can use observations of H^D" 1 ", in conjunction with Rd, to set limits on x e 
assuming the numbers of positive and negative charges are approximately equal (i.e., the core gas 
is neutral). Given t he abundances o f N2H" 1 ", N2D + and H2D + in B2, we can thus derive a lower 
limit on x e following iMiettinen et al.l (|2009l ): 
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x e > A(N 2 H+) + X(N 2 D+) + X(H+) + A(H 2 D+) 



(9) 



where we neglect multiply deuterated forms of H3 . As discussed in ^4.4.21 predictions of the o/p- 
H 2 D + ratio vary from ~ 0.1 to > 1, and we use a moderate o/p-H 2 D + = 0.5 to e s timat e X(H 2 D + ) 
from X(ortho-H 2 D + ). In a simple, steady state analytical model, Crapsi et al. ( 20041 ) relate Rd 



(where R D = iV(N 2 D+) / iV(N 2 H+)) to the relative H 2 D+-H+ abundance, r = [H 2 D+]/[H+], and 
also neglect mu l tiply deuterated forms of H3" such that Rd ~ r /(3 + 2r). In an updated result, 
Miettinen et ajj (|2009i ) find i?o « (r+2r 2 )/(3+2r+r 2 ). Using the lMiettinen et al.l relation we then 
find a mean r = 0.07 across Oph B2, and a lower limit x e > 1.5 x 10~ 8 through Equation^ We find 
small variations in the lower limit across B2, with a larger limit x e > 2.7 x 10~ 8 towards Elias 32 
and higher values in general near the 850 /im continuum peak. Detailed modeling, in conjunction 
with a measurement of the CO depletion factor in B2, is needed to improve this analysis of x e 
beyond a lower limit only. 

Low mass protostars emit x-rays, which could potentially increase the local ionization fraction 
and consequently decrease locally the deuterium fractionation of gas in Oph B2. In t heir study of the 
x-ray emission from young protostars in the central Ophiuchus region with ROSAT. ICasanova et al. 
(|1995l ) detected both Oph B2 protostars in the 1 — 2.4 keV range, finding luminosities (adjusted 
for extinction) L x = 30.7ergs _1 and L x = 29.2 ergs -1 for E32 and E33, respective ly. We can then 



estima te the x-ray ionization rate due to the protostars, Q x , at a distance r following lSilk fe Norman 
(|1983l . Equation 2). Prom Figure [T2l we find trends in Rd and X(ortho-H 2 D + ) to projected 
protostellar distances of order a few xO.Ol pc. Assuming a constant, uniform density n = 10 5 cm -3 , 
at r = 0.02 pc we find ( x ~ 9 x 10~ 17 s -1 , and at r = 0. 04 pc we find Cr ~ 8 x 10~ 18 s -1 . The 
fractional ionization based on these rates, again following ISilk &: Normanl . is x e ~ 9 x 10 -8 and 
i e ~3x 10~ 8 at r = 0.02 pc and r = 0.04 pc, consistent with our lower limits calculated above. 
Clumpy, rather than uniform gas could allow x-rays to penetrate further into the core, potentially 
increasing the fractional ionization t o greater distances. The st andard ionization rate due to cosmic 
rays alone Ccr = 1-3 x 10~ 17 s -1 (Her bst Sz Klempererl Il973l ). leading to a fractional ionization 
i e ~4x 10 -8 for n = 10 5 cm" " 3 (lMcKeelll989F " Our simple calculation shows that at distances 



comparable to where we find variation in the deuterium fractionation in Oph B2, the protostellar 
x-ray ionization rate and resulting fractional ionization may be comparable to or greater than 
expected due to cosmic rays alone. This result suggests that x-ray ionization is a viable mechanism 
to depress the deuterium fractionation close to embedded protostars. Further work is needed to 
directly measure the fractional ionization of dense gas near protostars, and to quantify the effect 
of a greater ionization rate on the chemical network of deuterated species in high density gas. 



5.3.4- Other possibilities 
While Figure [12] suggests that the distribution of enhanced deuterium fractionation of H^" and 
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N2H + in Oph B2 is tied to the presence of the Class I protostars associated with the core, other 
factors may influence the deuteration chemistry in this complex environment. For completeness, 
we briefly mention two other possibilities. First, line-of-sight inhomogeneities in Oph B2 would 
impact our estimate of the deuterium fractionation, as our analysis assumes constant excitation 
and abundance conditions along the line of sight. Detailed physical and chemical modeling may be 
useful in probing the effects of possible inhomogeneities, but is beyond the scope of this paper due 
to the complexity of Oph B2. Second, differences in the evolutionary timescale may be important in 
Oph B2, and in other multiple-star-forming objects. While several protostars have already formed 
out of the core gas (particularly in the south), there remain multiple dense, starless clumps within 
B2 which may yet form stars. The timescale on which different parts of the core condense and 
collapse may impact the dense gas chemistry and produce variations in the deuterium fractionation 
across Oph B2. 

In summary, there a number of different mechanisms by which the deuterium fractionation of 
N2H + and may be decreased near young protostars. We feel that a local increase in the gas 
temperature is the most likely explanation of this trend, but our simple calculation of the expected 
increase in the ionization fraction above the standard cosmic ray value is intriguing and worth 
further study. 

6. Summary 

In this paper we have presented the results of N 2 D+ 3-2, N 2 H+ 4-3 and ortho-H 2 D+ In - li 
mapping of the cluster-forming Ophiuchus B Core. In general, our results show that care must be 
taken when using deuterated species as a probe of the physical conditions of dense gas conditions 
in star-forming regions. We summarize in detail our results below. 

1. Significant emission from the deuterated species N 2 D+ and H 2 D + extends over several 
square arcminutes in B2. In particular, the extent of the H 2 D + emission is the largest yet mapped. 
Both species' integrated intensity distributions are offset from the continuum emission, and avoid 
entirely parts of B2 near embedded protostars and the peak of continuum emission from cold dust. 

2. Through fitting of the spectral lines, we find greater N 2 D + line widths near protostars 
embedded in B2, suggestive of protostellar influence on the dense gas. Complicated line structures 
are found, including blue line shoulders, in all three tracers. We estimate the non-thermal line 
widths for the species observed here and compare with NH3 and N 2 H + results from Paper I and 
Paper II. In B2, non-thermal line widths decrease with increasing critical density of the molecular 
tracer, but remain transonic or only slightly subsonic even at high densities. This result is in 
contrast with the nearly thermal line widths found in many isolated cores. 

3. We test a method to estimate the gas temperature by equating the non-thermal line widths 
of N 2 D + and H 2 D + and solving for Tk- We find both realistic (Tk ~ 11 K) and unrealistic 
(Tk < OK and Tk > 25 K) values in B2. Complicated line structure in B2 likely affects this 
analysis, but we expect good results can be found in regions with simple velocity structure. 
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4. The deuterium fraction of N 2 H+ R D = iV(N 2 D + )/N (N 2 H+) ~ 0.03 on average where 
N 2 D + is detected. Small-scale features of enhanced Rd are seen in B2, but most do not correlate 
with previously identified continuum clumps or N 2 H + clumps. This Rd is consistent with previous 
results in protostellar cores, and is also within the lower range of values found for starless cores. 

5. The average ortho-H 2 D+ abundance, A"(ortho-H 2 D + ) ~ 3 x 10~ 10 , given an assumed 
excitation temperature of 12 K. An accurate measure of the ortho- to para-H 2 D + ratio is needed 
to convert to total H 2 D + abundance, and has been found in previous studies to vary between 0.1 
and > 1. Future observations of para-H 2 D + transitions, if possible, are needed to constrain this 
value in dense cores. 

6. An anti-correlation is found with Rd and X(ortho-H 2 D + ) and projected distance to the 
nearest embedded protostar in B2. This relation indicates that the embedded protostars in B2 
are affecting the internal chemistry. Possible mechanisms are through gas heating, evaporation 
of CO off dust grains, CO mixing caused by the B2 outflow, or by increasing the local ionization 
fraction. Gas temperatures determined through NH3 analysis in Paper I are < 18 — 20 K required to 
evaporate CO and reduce deuterium fractionation, but NH3 may not be a good tracer of conditions 
in the high density gas in Oph B2. We calculate that the x-ray flux of the embedded protostars in 
B2 could increase the fractional ionization of the core gas beyond that expected solely due to cosmic 
ray ionization to distances ~ 0.02 — 0.03 pc, similar to the distances over which we find variation 
in the deuterium fractionation. We suggest that a temperature increase is the dominant parameter 
in this trend, but further research on the effect of young stars and their outflows on both the local 
ionization fraction and the chemistry of deuterated species in protostellar cores is needed. 

7. We use the observed Rd and X(ortho-H 2 D + ) to estimate a lower limit on the ionization 
fraction in B2, and find x(e) > 1.5 x 10~ 8 given an assumed ortho- to para-H 2 D + ratio of 0.5. 
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Table 1. Observed Species, Transitions and Rest Frequencies 



Species 


Transition 


F{ 


F' 


Fi 


F 


Frequency 
GHz 


N 2 D+ 


3-2 


4 


5 


3 


4 


231.321906 a 


N 2 H+ 


4-3 


5 


5 


4 


4 


372.672494 b 


H 2 D+ 


In — lio 










372.421385 c 



Gerin et all (j200lf ) 



Pagani et all (|2009j ) 



Amano & Hiraol (|2005l ) 
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a (2000.0) 



Fig. 1. — a) Integrated N 2 D + 3-2 emission towards Oph B2 (greyscale) in Tg units. Data have 
been convolved to a final 18" beam (FWHM). Black contours show 850 p continuum emission 
in increments o f 0.1 J y beam -1 . Squares show the locations of continuum clumps identified by 
Jprgensen et al.l (|2008l ). The continuum clump B2-MM8 is identified. In both plots, stars indicate 
the positions of Class I protostars (with Elias 32 and 33 identified) and the 18" beam is shown at 
lower left, b) Integrated N 2 H + 1-0 emission towards Oph B2 observed with the Nobeyama 45 m 
Telescope at 18" (FWHM) resolution (greyscale). Black and grey contours begin at 3Kkms _1 
(Tmb) and increase by lKkms -1 . Triangles show the peak locations of N 2 H + clumps identified 
with CLUMPFIND, and the continuum clump B2-MM8 is shown by the square. White contours trace 
N 2 D + emission from a), beginning at 0.45 K km s -1 (Tg) and increasing by 0.15 K km s -1 . 
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Table 2. Mean NH 3 , N 2 H+, N 2 D+ and H 2 D+ v LSR and Av in Oph B2 



Species Transition 



VLSR 
kms -1 



Av 
kms -1 



cuT 









Mean 


RMS 


Min 


Max 


Mean 


RMS 


Min 


Max 




NH 3 


(J,K) = 


= (1,1) 


4.14 


0.19 


3.85 


4.58 


0.87 


0.16 


0.42 


1.11 


2 x 10 4 a 


N 2 H+ 


J=l- 





4.04 


0.20 


3.61 


4.43 


0.71 


0.22 


0.38 


1.42 


1.4 x 10 5b 


H 2 D+ 


JK-xKx 
























= 111 - 


lio 


4.01 


0.24 


3.65 


4.50 


0.74 


0.19 


0.51 


1.36 


~ 10 5c 


N 2 D+ 


J = 3- 


2 


4.02 


0.19 


3.67 


4.44 


0.53 


0.20 


0.25 


1.36 


8 x 10 5 d 


N 2 H+ 


J = 4- 


3 e 


4.05 


0.22 


3.39 


4.43 


0.69 


0.18 


0.47 


1.14 


7.7 x 10 6b 



Note. — The velocity resolution of the NH3 observations was 0.3 kms . All other species were 
observed with 0.1 kms -1 velocity resolution or better. 



HCaselli et all (2002a|) 



nSchoier et al.l (n cr calculated using Einst ein A coefficients and collision rates j u i at 10 K from 
the Leiden Atomic and Molecular Database; 



2005) 



Caselli et all faOO® ) 



d Cal culated using Einst ein A coefficients and collision rates for N 2 H + 3-2 at 10 K ( Schoier et al 



2005 


)• 


Daniel et al. 


(2007^ 



coefficients and collision rates, critical densities of N 2 H + and N 2 D + are similar for transitions other 
than J = 1 — 0. 



e Values calculated from a Gaussian fit due to a lack of visible hyperfine structure. The reported 
line widths may be greater by ~ 10 — 25 % than the true values; see §4. H 
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Fig. 2.— H 2 D+ In - lio spectra observed with HARP at the JCMT towards Oph B2 with 15" 
FWHM. The grey contours show 850 /urn continuum emission observed at the same resolution at 
O.lJybeam -1 , 0.3Jybeam _1 and 0.5Jybeam _1 . Stars represent Class I protostars, with solid 
lines indicating the location if within a spectrum box. The velocity scale on the subplots runs 
from -lkms _1 to 9kms -1 , while the amplitude ranges from -0.25 K to 0.85 K {Tmb-, assuming 

1MB = 0.6). 
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Fig. 3.— N 2 H+ 4-3 spectra observed with HARP at the JCMT towards Oph B2 with 15" FWHM. 
The grey contours, stars and velocity scale on the subplots are as in Figure [2 The amplitude ranges 
from -0.3 K to 2.2 K (Tmb, assuming tjmb = 0.6). Four locations of interesting spectra, labeled A, 
B, C and D, are highlighted by dark outlines. 
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v lsr (km s 1 ) 

Fig. 4. — N2H + 4-3 spectra towards four locations in Oph B2, labeled A, B, C and D in Figure [3j 
30" north of E33 (A), towards E33 (B), 21" north-east of E33 (C), and ~ 42" north-east of E33 
(D). Overlaid on the spectra are fits to the full hyperfine structure of the N2H + 4-3 line (A, B and 
D) and a two-component Gaussian fit (C). 
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Fig. 5. — a) N2D + 3-2 line velocity or vlsr in Oph B2. Colour scale is in kms -1 . In both plots, 
only results from pixels where the S/N of the peak line intensity is > 5 are shown, indicated by 
the solid white contour. Contours show 850 /jm continuum emission at 15" (FWHM) resolution 
observed with the JCMT, in increments of O.lJybeam -1 . The 18" (FWHM) beam is shown at 
lower left. Stars indicate the positions of Class I protostars. The locations of the NH3 clump B2-A7 
and 850 ^m continuum clump 16275-2422 are also shown, b) Fitted Av in Oph B2, contours as in 
a). Colour scale is in kms -1 . 



-36- 



26 00 

^ 30 

o 

d 

g 27 00 

30 

28 00 

-24 28 30 

26 00 

^ 30 

o 

o 

g 27 00 
cv 

30 
28 00 

-24 28 30 



i 




y 



II ii ISS7B8I ill 



16 27 36 



33 



30 27 24 
a (2000.0) 



21 



16 27 36 33 30 27 24 
a (2000.0) 



21 



Fig. 6. — a, c) H2D + In — lio and N2H + 4-3 line velocity f LSR in Oph B2, respectively. Colour 
scale is in kms -1 . We show only results from pixels where the S/N of the peak line intensity is 
> 5. Grey contours show 850 //m continuum emission as in Figure [5j Stars indicate positions 
of embedded Class I protostars. b, d) H2D + and N2H + line width Av in Oph B2, respectively. 
Colour scale is in kms -1 . 
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Fig. 7.— a) Total N 2 D+ column density iV(N 2 D+) (cm 2 ) in Oph B2 determined using Equation 
[3] (greyscale). Values shown have been divided by 10 11 . In both plots, only results from pixels 
where the S/N of the peak line intensity is > 5 are shown, indicated by the thick black contour. 
Grey contours show 850 /um continuum emission as in Figure [5j The 18" (FWHM) beam is shown 
at lower left. Stars indicate the positions of Class I protostars. b) Fractional abundance of N2D + 
relative to N2H + , Ftp, in Oph B2. The scale has been truncated to show the variation in central 
B2; the maximum Ftp = 0.16 towards the eastern tip. 
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Fig. 8. — Impact of the assumed excitation temperature, T ex , on returned H2D + In — lio line 
opacity (top), partition function Q rot (middle), and column density, iV(ortho-H2D + ), given a line 
brightness temperature Tmb and line width Av similar to observed mean values in Oph B2. At 
low T ex , t and iV(ortho-H2D + ) vary strongly with T ex , but the effect of changes in T ex at higher 
temperatures becomes small. 
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Table 3. Impact of T ex on r, Q ro t and -/V(ortho-H2D + ) 





r 


Qrot 


7V(ortno-H 2 D + ) 


K 






10 J cm z 


7 


0.41 


9.7 


1.60 


8 


0.27 


10.0 


1.11 


9 


0.20 


10.2 


0.86 


10 


0.15 


10.5 


0.71 


11 


0.12 


10.8 


0.61 


12 


0.10 


11.0 


0.54 


13 


0.09 


11.3 


0.48 


14 


0.08 


11.5 


0.45 


15 


0.07 


11.7 


0.42 



Note. - Values were calculated 
for a line strength Tmb = 0.5 K and 
width Aii = 0.7 km s -1 . 



5 10 15 012345 




a (2000.0) a (2000.0) 

Fig. 9. — a) Ortho-H2D + column density in cm~ 2 (greyscale). A r (ortho-H2D + ) values have been 
divided by 10 12 . Grey contours show 850 /im continuum emission as in Figure EJ b) Fractional 
ortho-H2D + abundance X(ortho-H2D + ) (greyscale). Values have been divided by 10~ 10 . 
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Fig. 10. — a) Tk determined using Equation [7] for H^D" 1 " and N2D + emission. Note that a negative 
Tk was found for the blank pixel towards the submillimetre continuum emission peak because the 
observed H2D + Av > N2D + Av. Grey contours show 850 /im continuum emission as in Figure [5j 
b) Tk determined using Equation [7] for H2D + and N2H + 4-3 emission. 
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Fig. 11. — Top left: H2D + In — lio spectrum towards B2-MM8 overlaid with a Gaussian fit. 
Top right: N2D + 3-2 spectrum towards B2-MM8 overlaid with a hyperfine structure fit (note 
no individual hyperfine components are visible due to their small separation in velocity and the 
relatively large line width) . Bottom: Returned H2D + line widths Av versus N2D" 1 " Av. Each point 
represents a 15" pixel. An equality relation is shown by the dotted line. Despite appearing wider, 
in most locations the N2D + Av is narrower than that of H2D + due to the presence of overlapping 
hyperfine lines. The two points where N2D + At; > H2D + At; are located at and immediately 
adjacent to the 850 /um continuum peak and B2-MM8. At these locations, temperatures calculated 
via Equation [7] become negative and unphysical. Uncertainties shown are returned from the line 
fitting routines. Uncertainties in the line widths are shown. 
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Fig. 12. — Variation in Rp, A(N2H + ), X(ortho-H2D + ) and Tk with H2 column density and the 
projected distance to the nearest embedded protostar in Oph B2. In each plot, points shown 
represent pixels matched to the FWHM beam of the respective datasets (15" for Tk and H2D + , 
and 18" for N2H + and N2D+). Typical uncertainties are shown for a single pixel in each plot. 
An uncertainty of ~ 20 — 30% in the A. (ortho-H^D -1- ) values is not shown. Only pixels with S/N 
> 5 are plotted. In (a), no trend is found in Ru versus logiV(H2). The Rd data are, however, 
consistent with an increasing deuterium fractionation with distance from an embedded source (b), 
albeit with significant scatter. As discussed in Paper II, a strong correlation is seen in X(N2H + ) 
versus logiV(H2) (c), but no trend is seen with respect to projected protostellar distance (d). 
A(ortho-H2D + ) shows no trend with log A(H2) (e). With respect to projected protostellar distance, 
the X(ortho-H2D + ) plot shows the data do not scatter about a central point (f), but instead show 
a trend of increasing A(H2D + ) to d ~ 0.04 pc, and a smaller, constant H2D + abundance at further 
distances. In both (e) and (f), the points identified with squares represent pixels in the northeast 
portion of the B core. A linear fit to the data, omitting the highlighted pixels, is shown by the 
dotted line. A linear fit to the data is shown by the dashed line. No clear correlation is found in 
Tk versus logiV(H2) (g) or projected protostellar distance (h). 



